Transfer ionization and its sensitivity to the ground-state wave function 
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We present kinematically complete theoretical calculations and experiments for transfer ionization 
in H + +He collisions at 630 keV/u. Experiment and theory are compared on the most detailed level 
of fully differential cross sections in the momentum space. This allows us to unambiguously identify 
contributions from the shake-off and two-step-2 mechanisms of the reaction, ft is shown that the 
simultaneous electron transfer and ionization is highly sensitive to the quality of a trial initial-state 
wave function. 



INTRODUCTION 



Double ionization of Helium is the benchmark system 
to study electron-electron correlation in many electron 
systems and to test state-of-the-art theories. One [HQ, 
two [3] and multiphoton [i| double ionization has been 
explored in great detail experimentally and theoretically. 
For single photon ionization it has been shown that only 
fully differential cross sections reveal the mechanisms for 
double ionization as shake-off (SO), the so called two-step 
processes (TS) Q and the quasi free mechanism (QFM) 
3, 0]. For ion impact the state-of-the-art is much less 
satisfactory. Even for single ionization by fast particle 
impact unresolved discrepancies between theory and ex- 
periment remain @~Q ■ This is even more problematic for 
transfer ionization (e. g. H + + He — > H° + He 2+ + e~) 
where projectile bound states have to be included 13- 13] . 
So far no quantum theory is available which is able to re- 
produce the full 3 dimensional momentum distributions. 
A theory which is capable to predict all the experimen- 
tal observables for transfer ionization (TI) is a particular 
worthy goal as there are many indications in the litera- 
ture that transfer ionization is an extremely interesting 
channel, which rich features that cannot be accessed by 
photon, ion or electron impact double ionization or in 
strong laser pulses. 

A physical explanation how a target electron can be 
captured into a bound state of a fast moving projec- 
tile within the single-interaction scenario was given by 
Oppcnhcimer, Brinkmann and Kramers (OBK) 171 Il8|. 
In the OBK approximation the electron transfer pro- 
ceeds via a momentum space overlap of the initial target 
and final projectile wave function, which are displaced 
by the projectile velocity v p . Thus kinematical capture 
at velocities above the Bohr velocity relies on the high- 
momentum components in the ground state wave func- 



tion. Therefore the kinematical capture steeply decreases 
with projectile velocity (a cx v~ 12 ) [l7l. Hij. 

Emission of the second electron can take place via 
shake-off (SO) due the sudden removal of its correla- 
tion partner in the bound state [3, EoT - E^ . This suggests 
to treat the shake-off following a kinematical capture 
in analogy to shake-off following photoionization 23| . 



since in the OBK approximation as well as in the pho- 
toionization process the first electron is removed from 
a certain velocity component of the ground state (see 
13, E3, E3, 0, III] ) . Experiments of Mergel et al. [13, EI | 
raised the question whether higher angular momentum 
components in the ground-state wave function (so-called 
non Is 2 contributions), though contributing only about 
2 % of the total wave function, play a significant role in 
producing the observed momentum distributions in the 
continuum. This has been later supported by calcula- 



tions of Godunov and coworkers [20|, [21|, |24|, |26j . Also a 
second interaction projectile-electron, e. g. binary en- 
counter (BE), can take place, leading to emission of the 
second electron. This process is usually termed indepen- 
dent two-step 2 (TS2) mechanism 27 1. 

Recently Voitkiv predicted a further new higher order 
mechanism for transfer ionization. This process is sim- 
ilar to radiative capture [28| . but here the energy from 
the relative motion between projectile and target frame is 
transferred to the second target electron, which is emit- 
ted opposite to the incoming projectile 13, 13]. 

At high velocities in addition the classical Thomas pro 
cess, which includes electron-electron interaction 
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22]- [3l|, contributes. It is of the second order, which be- 



comes significant only at much higher projectile energies 
and scattering angles and is not in the focus of this paper. 

In this letter we present for the first time the exper- 
imental distribution of the momentum of the escaped 
electron in the scattering plane and the corresponding 
calculations in the plane wave first Born approximation 
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(PWFBA) on the level of fully differential cross sections. 
This is the most sensitive test of theory possible. The 
present calculations yield an unprecedented insight into 
the physical mechanisms at play as they allow to change 
the initial state and selectively switch on and off the dis- 
tinct mechanisms discussed above. The advantage of 
the present calculations is that the initial state can be 
changed and the different mechanisms can be selectively 
switched on and off. By comparing these calculations 
to our high-resolution experimental data, we separate 
ionization due to shake-off (A\ + A3) or binary collision 
(A 2 ) (see Fig. 1) leading to distinct islands in momentum 
space. One of the results is that these data are extremely 
sensitive to the initial-state correlation. 
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turn of the projectile originates from the nucleus-nucleus 
respectively electron-proton scattering, the only source of 
transverse momentum in A 2 originates from the binary 
collision between the proton and the emitted electron. 
Therefore the contributions from A2 will concentrate at 
unrealistically small scattering angles. 

The captured electron always moves in the direction 
of the initial proton. In the case of strong correlation of 
electrons in the helium target, the electron ejected via 
shake-off should move in the opposite direction, while 
the mechanism A 2 is responsible predominantly for the 
ejection of this electron in the forward direction. It was 
shown in j32j that the term A 2 can give a rather big and 
even the leading contribution to the differential cross sec- 
tion in the case of TI processes. Therefore we can expect 
to observe noticeable distributions both in forward and 
backward directions when the trial helium ground-state 
wave function is well correlated and contains higher an- 
gular momenta (non Is 2 part). 
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FIG. 1: Schematic presentation of the A\ + A3 and A2 terms. 
Ai (OBK) describes the collision between electron and proton 
followed by the capture of this electron by the projectile. The 
second electron is released due to rearrangement in the he- 
lium, known as shake-off (SO). The A3 amplitude represents 
first an interaction between the projectile and the helium nu- 
cleus followed by electron capture. Similarly to Ai the second 
electron is also released due to the sudden rearrangement in 
the helium. Term A2 describes the classically termed TS2 
amplitude. First the proton knocks-off a target electron into 
the continuum, followed by capturing the remaining electron 
from the helium. 



In detail (see Fig. 1) the term A\ (OBK) describes 
the collision between the proton and a target electron, 
which then is captured by the projectile. The second 
electron is released via SO. The collision between the 
proton and the target electron deflects the proton, lead- 
ing to small scattering angles <0.5 mrad. Also in A 3 
the electron emission takes place due SO, following the 
electron capture. In contrast to A\ here, the capture 
proceeds after a nucleus-nucleus scattering of target and 
projectile followed by the electron capture. This term 
accesses also the larger scattering angles and dominates 
above 0.5 mrad. 

In term A 2 the first interaction directly knocks-off a 
target electron and the second electron is captured. Even 
commonly termed two-step 2 (TS2) mechanism [27[, this 
process can be described in a first Born approximation, 
as strictly speaking the capture can be of zero's order. 
Contrary to terms A\ + A3 where the transverse momen- 



EXPERIMENT 

We have used the COLTRIMS technique [H-HH to de- 
termine the momentum vectors of all final-state products. 
The experiment was performed at the Van de Graaff ac- 
celerator of the Institut fur Kcrnphysik at the University 
of Frankfurt. The projectile beam (H + ) was collimated 
to a size of about 0.5 x 0.5 mm 2 at the target. 15 cm 
upstream of the target, a set of electrostatic deflector 
plates cleaned the primary beam from charge state im- 
purities. The proton beam intersects with a supersonic 
Helium gas jet (density of 5 x 10 11 atoms/cm 2 and a di- 
ameter of 1 mm). About 15 cm downstream a second set 
of electrostatic deflector plates separate the final charge 
state, thus only the neutral projectiles (H°) hit a posi- 
tion and time sensitive multichannel plate (MCP) detec- 



tor |36|, |37| yielding the projectile deflection angle and 
the time zero of the collision. The recoil ions were accel- 
erated by a weak electrostatic field of 4.8 V/cm in the in- 
teraction region and detected on a 80 mm MCP-detector 
with delay-line anode. To optimize the resolution, a three 
dimensional time and space focussing geometry 38, 39| 
was used for the spectrometer. A momentum resolution 
of 0.1 a. u. was achieved in all three directions. The elec- 
trons were guided by a magnetic field of 15 Gauss and 
accelerate d by the same electric field in a time focussing 
geometry [4QJ onto a multi channel plate detector of 120 
mm active diameter. A three-particle coincidence (H°- 
He 2+ -e) was applied to record the data event by event. 
From the positions of impact on the detectors and the 
time-of-flight we can derive the initial momentum vectors 
of the He ++ and the electron. Energy conservation was 
used for offline background suppression. Furthermore the 
high resolution data allowed to distinguish data where 
the neutral projectile (H°) is found in an excited state 
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from those, where the hydrogen is in its ground-state. 
Only these latter ones are presented in this letter. 



THEORY 

In this study we examine very small scattering angles 
of the outgoing hydrogen (0 < 9 p < 0.5 mrad). It leads 
to a practically zero recoil-ion velocity in the laboratory 
frame during the process, so we treat it as being at rest. 

In the presented experiments, the scattering plane 
formed by the momentum vectors p p (z-axis) and pn is 
fixed in space, and we choose its polar angle to be = 0. 
The Double Differential Cross Section (DDCS) takes the 
form 



2m 2 
(2tt) £ 



J 8pd8 p J dk y \A, +A 2 + A 3 \ 2 . (1) 



Here m is the proton mass, the scattering angle 

domain and (k x , k y , k z ) the electron momentum compo- 
nents. We calculate the DDCS depending on (k z ,k x ) 
electron momentum distribution in the scattering plane. 

In theoretical calculations we use two trial ground- 
state helium wave functions. One is the loosely corre- 
lated Is 2 Roothaan-Hartree-Fock (RHF) wave function 
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no angular correlation) with a rather poor ground- 
state energy of -2.861680 a. u. Another one is the highly 
correlated wave function given in 42| with a ground-state 
energy of -2.903721 a. u. being very close to the experi- 
mental value. 



RESULTS AND DISCUSSION 

In Fig. 2 we present experimental electron momen- 
tum distributions and theoretical results in the scatter- 
ing plane (the x component of the vector pn is positive 
here). In the left column (Fig. 2a,c,e) only events for a 
small projectile scattering angle 9 P <0.25 mrad are se- 
lected and in the right column (Fig. 2b,d,f) events for 
0.25< 6 p <0.45 mrad. The experimental data in Fig. 2a 
show that at these small scattering angles, the electron 
is predominantly emitted in backward direction, while at 
larger scattering angles, corresponding to more violent 
collisions, shown in Fig. 2b, the electron is found close 
to the binary encounter ridge indicated by the red dashed 
line. 

The results using a Is 2 trial helium wave function 
are shown in Fig. 2c, d. Particularly for small 9 p the 
momentum distribution is very similar to the binary- 
and recoil-peak structures (forward and backward) well 
known from electron impact ionization experiments |43| ■ 
For comparison, the separation of the individual con- 
tributions/processes is shown in Fig. 3; the expected 
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FIG. 2: Experimental and theoretical data for 630 keV 
H+/He collisions for 6 P <0.25 (left column) and 0.25< 
9 P <0.45 mrad (right column). (a,b) are the experimental 
data and the red dashed line represents the binary encounter 
ridge. (c,d) calculations using a helium Is 2 trial wave func- 
tion, while (e,f) uses a highly correlated helium wave function 
with angular momentum, all including terms A 1 + A2 + A 3 



electron momentum distribution for the shake-off-process 
(Ai + ^3-term) in the case of loosely correlated helium 
wave function is shown in the top row (Fig. 3a, b), while 
in Fig. 3c, d (lower row) only the sequential TS2 mech- 
anism (A 2 -term) is taken into account. The shake-off 
exhibits a perfectly isotropic behavior, as expected for a 
ls 2 -state with zero angular momentum. This term has 
however a visible influence on the coherent sum of con- 
tributions Al, A2 and A3 (Fig. 2c, d) in spite of the 
dominance of the A2 (TS2) term in whole. It changes 
the binary/recoil peak ratio, while conserving the general 
features of forward and backward contributions leaving 
the general distributions to be similar. Comparing the 
experimental and theoretical results (couples in Fig. 2a, c 
and Fig. 2b, d), we find notable differences. 

The agreement improves considerably for a well- 
correlated helium wave function with radial and angular 
ee-correlations. In Fig. 2e,f the results of our calculations 
are shown for small and intermediate 8 P values. And they 
are split into the different contributions in Fig. 4. It 
can clearly be seen that the shake-off terms in Fig 4a, b 
show an asymmetric emission pattern, peaking in back- 
ward direction. A binary /recoil peaklike structure is also 
clearly visible for the A2-term using a correlated wave 
function (Fig. 4c, d). The coherent sum (Fig. 2e) also 
exhibits two clearly distinct peaks pointing forward and 
backward along the z-axis. This structure is considerably 
rotated clockwise compared to the one shown in Fig. 2c. 



FIG. 3: Calculation for a helium Is 2 trial wave function sep- 
arated in SO (a,b) and TS2 (c,d) for 6 P <0.25 (left column) 
and 0.25< 8 P <0.45 mrad (right column). 



FIG. 4: Calculation for a highly correlated helium wave func- 
tion with angular momentum separated in SO (a,b) and TS2 
(c,d) for 6 P <0.25 (left column) and 0.25< 8 P <0.45 mrad 
(right column). 



While the calculation in Fig. 2c suggest a more equal 
contribution of A\ + A3 vs. Ai the experiment shows 
predominantly the backward contribution from A\ + A3. 
This is clearly a consequence from the ovcrestimation of 
the Ai term at small scattering angles caused by the 
missing nuclear-nuclear deflection term in A2 mentioned 
above. A calculation in the second Born approximation, 
which would account for nucleus-nucleus interaction and 
the momentum kick from the electron transfer would def- 
initely solve this. It has to be noted that the process 
proposed by Voitkiv [29[ also predicts electrons emitted 
in backward direction in agreement with experiment [3fj| . 
So far no fully differential cross sections have been pub- 
lished for the mechanism predicted in 0| . We note that 
this mechanism is not implicitly included in our PWFBA 
calculations 44j . Our present work suggests that such an 
additional mechanism beyond PWFBA is not necessary 
in order to explain the backward electron emission in TI 
at the present collision velocity. 

Comparing now the experimental and theoretical re- 
sults (couples in Fig. 2a, e and Fig. 2b, f), we find good 
agreement, which is yet less satisfactory for at larger pro- 
jectile scattering angles 0.25< 8 P <0.45 mrad (right col- 
umn). At larger scattering angles we expect also the 
second Born term to contribute which is not included in 
our present TI calculations. 



son of a loosely and a strongly correlated wave function 
for the initial state confirms the high sensitivity of the ex- 
periment to the subtle features of the initial state wave 
function. Even better agreement especially for the for- 
ward emitted electrons, TS2 contribution (A^-term), can 
be expected for calculations in the second order. 
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CONCLUSIONS 

In conclusion, wc presented highly differential theory 
(PWFBA) and experimental data from a kinematical 
complete experiment on transfer ionization in proton- 
Helium-collision at 630 keV/u. The observed split- 
ting into forward and backward emission originates from 
two different contributions, the A2-term (TS2, electron 
knock-off) and the A\ + A3-term (shake-off). Compari- 
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